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Abstract The lifetime of solar-like stars, the envelope 
structure of more massive stars, and stellar acoustic fre- 
quencies largely depend on the radiative properties of 
the stellar plasma. Up to now, these complex quan- 
tities have been estimated only theoretically. The de- 
velopment of the powerful tools of helio- and astero- 
seismology has made it possible to gain insights on the 
interiors of stars. Consequently, increased emphasis is 
now placed on knowledge of the monochromatic opacity 
coefficients. Here we review how these radiative prop- 
erties play a role, and where they are most important. 
We then concentrate specifically on the envelopes of 
/3 Cephei variable stars. We discuss the dispersion of 
eight different theoretical estimates of the monochro- 
matic opacity spectrum and the challenges we need to 
face to check these calculations experimentally. 
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After more than two decades of study, helioseis- 
mology has demonstrated the capability to probe 
the solar interior thanks to the seismic instruments 
of the SOHO satellite and ground-based networks 
( Christensen-Dalsgaard 2004[ Guzik, Watson & Cox 
2005 ^Basu and Antia |2008[ |Turck-Chieze fc Couvidat 
,2010ar 



The new space missions (SDO, PICARD) are ori- 
ented toward a better understanding of the Sun's in- 
ternal magnetism and its impact on the Earth's atmo- 
sphere. In this context, the existence of a fossil field in 
the Sun's radiative zone is a missing bridge that mo- 
tivate us to better estimate the competition between 



microscopic and macroscopic physics ( 


Duez, Mathis & 


Turck-Chieze ||2010 


Turck-Chieze et al.||2010b 


). 



In addition to the international solar research com- 
munity, more than 350 scientists from all over the 
world are working on the data from the COROT 
and Kepler spacecraft dedicated to exoplanets and 
astero-seismology (see http://smsc.cnes.fr/COROT for 
COROT and http: / /astro.phys.au.dk/KASC / for Ke- 
pler). The development of these disciplines justifies en- 
suring the correctness of the many fundamental proper- 
ties of the plasma (detailed abundances, reaction rates, 
opacity coefficients, etc.). 

The new High Energy Density laser facilities and 
the Z pinch machine now enable access to the radiative 
properties of plasmas of stellar type in the laboratory 



during pico and nano second timescales; see (Bailey et 
al. ||2007[ [Bailey et al.|[2009{ [Turck-Chieze et al.||2009[ ), 

and references therein. This review is divided into three 
parts: Section 2 is dedicated to solar-type stars. Sec- 
tion 3 to the more massive stars, and Section 4 to the 
study of the /? Cephei variable star case and a com- 
parison between eight calculated iron opacity spectra. 
The challenges faced by the experiments are described 
in Section 5. 
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2 Radiative properties of solar-type stars 

Stars produce both nuclear energy through nuclear re- 
actions, and thermal energy due to the internal temper- 
ature that exceeds everywhere some thousand degrees. 
The central plasma is heated to more than a dozen mil- 
lion K, and the heat is transported by convection or 
radiation phenomena. 

The low-mass stars (< 1.5Mq) have an internal ra- 
diative zone which disappears in the region where the 
photons interact so strongly with the partially ionized 
constituents that convection becomes the most efficient 



highlight the important effect of the heavy elements on 
the sound speed profile ( Turck-Chieze et al.|2004| Bah- 



way to transport the energy; see Figure 32 of Turck-| 



Chieze et al. (1993) or Turck-Chieze et al. (2010b I. 



In fact, this boundary between radiative transport 
and convective transport differs from star to star de- 
pending on the star's mass, age, and composition. Hy- 
drogen and helium are completely ionized in these ra- 
diative zones due to their high temperature and den- 
sity, but in the central conditions, the heaviest elements 
(iron and above) are generally partially ionized and 
their contribution to the opacity coefficients is rather 
high, 20% in the solar case; see equation 3 and figures 
in Turck-Chieze et al. (2009). Nevertheless, their mass 
fraction is of the order of 10"'^. In the rest of the ra- 
diative zone, since both the temperature and density 
decrease with radius, by typically a factor 15 and 100 
respectively, most of the heavy elements down to oxy- 
gen become partially ionized, so their relative contri- 
bution increases the Rosseland mean opacity kh by a 
factor up to 100 at the base of the convection zone due 
to the multiple photon interactions on the electronic 
bound layers of more and more elements. 

As a direct consequence, the lifetime of a star of a 
given mass strongly depends on the amount of heavier 
elements (carbon to uranium) present in the star. The 
total element mass fraction is generally divided into X, 
Y, and Z, with X and Y corresponding to the mass 
fraction of hydrogen and helium, respectively, and Z 
corresponding to the remainder. The lifetime of a 0.8 
Mq star with Z=0.001 is about 14 Gyr, while the hfe- 
time of a 0.8 Mq with Z=0.02 is about 22 Gyr. When 
the heavy element contribution is small, the star con- 
tracts more quickly, the central temperature is higher, 
the opacity coefficients are different, and the nuclear 
reactions become more efficient. 

This example shows the importance of verifying some 
detailed contributions to the total Rosseland mean 
opacity values to properly estimate how kr changes 
from star to star depending on its mass and on its de- 
tailed composition. 

The best known case is the Sun where where our 
efforts to incorporate new updates of the composition 



call et al.|[2005'; 'Asplund et al."2009^ 'Turc k-Chieze fc 



Couvidat 



2010a: Turck-Chieze et al. 
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2010b|. The dif- 
ference between the standard model sound-speed pre- 
dictions and the observed sound speed is larger than 
the error bars coming from the knowledge of the com- 
position and from the known difference between OPAL 



and OP mean Rosseland opacity tables (Badnell et al 



20051 



While for the solar conditions the Rosseland mean 
opacities agree in general among different calculations, 
the detailed comparison, element by element, shows 
rather large differences, as mentioned by Blancard et al. 
in Turck-Chieze et al. 2010c. So it is useful to compare 
opacity calculations of a specific element to laboratory 
experiments; then the extrapolation to other stars of 
different composition will be more reliable. The first 
experiments have been performed at the Sandia Z ma- 
chine to check the iron contribution at a temperature 
and density not far from those at the base of the solar 
convection zone (Bailey et al. 2007 Bailey et al.|2009 ). 
Extensive measurements on the NIF or LMJ are needed 
for several well-chosen elements and temperatures be- 
tween 193 and 400 eV to validate opacity calculations 
in the radiative zones of solar-like stars. 

Indeed the region just around the transition is clearly 
the most difficult region due to the interplay of micro- 
scopic and macroscopic phenomena (turbulence, rota- 
tion and magnetic fields). In fact, in contrast with the 
rest of the radiative zone, the boundary of the internal 
radiative zone depends directly on the detailed abun- 
dances of oxygen, iron, and neon and on the associ- 
ated Rosseland mean opacity coefficients. It depends 
also indirectly on the turbulence of these layers associ- 
ated with the fact that the stars are generally rotating. 
So the thermal instability produces a strong horizon- 
tal shear between a differentially latitudinal rotating 
convective zone and a reasonably flat rotation in the 



radiative zone (Turck-Chieze et al. 2010b). Moreover, 



the magnetic field in these layers probably plays a cru- 
cial role in the solar dynamo at the beginning of the 
11-year cycle. 

In the sub-surface layers of these stars, radiation and 
convection are in competition and the heavy elements 
are less important because most of the elements are 
neutral or even in molecules. In the corona, the plasma 
is heated to a million degrees, so the heavy elements 
play a crucial role, but the very low density plasma is 
no longer in LTE (local thermodynamic equilibrium), 
so non-LTE calculations are needed to explain the ob- 
served lines ( Colgan et al.||2008 ). 
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stellar structure 




Fig. 1 a: Temperature gradient ratio between a radiative 
or convective region versus temperature for a 10 Mq star, 
b: Temperature dependence of the total Rosseland mean 
opacity for OP calculation (solid blue line) and OPAL cal- 
culation (dashed red line), c: Relative contribution of the 
heavy elements to the total opacity, d: Relative contribu- 
tion of the light elements to the total opacity. 

3 Radiative properties of more massive stars 

When the stellar mass increases, the central tempera- 
ture increases and consequently carbon, oxygen, nitro- 
gen can penetrate the Coulomb barrier. In these con- 
ditions the hydrogen burning is no longer dominated 
by the weak interaction but by the strong interaction 
where these elements act as catalysts. The nuclear 
burning becomes more rapid, and the lifetime of these 
stars decreases from billions to millions of years, the 
central temperature increases strongly and the radia- 
tive diffusion has no time to take place. The core of 
these stars (> l.SM©) becomes convective, while the 
envelope is radiative. In the envelope, many elements 
are partially ionized, so the detailed calculation of the 
Rosseland mean opacity is even more complex because 
many electronic shells are open and the bound-bound 
processes are extremely important. 

Figure 1 illustrates the case of a 10 M0 star. All the 
species including iron are totally ionized in the center 
of this star due to high temperature, 34 MK, and a rel- 
atively low density of 8 g/cm'^. The opacity is not as 
high (Fig lb) as in lighter stars, but the luminosity due 
to the CNO cycle is so high that Wrad is higher than 
Vad (Fig la). When the ratio of gradients is greater 
than one, the region is convective; in the opposite case, 
the region is radiative. The rest of the star is radiative 
except in two small regions which correspond to the 
partially ionized contribution of helium and hydrogen 
at low temperatures (see Fig Id) and of the iron-group 
elements above 100,000 K (Fig Ic). Figure 2 shows that 
the differences between OP and OPAL are particularly 
important for Ni, Mn and Cr and that the contribu- 
tion of Cr and Ni could be very important in the case 
of OP. Therefore, we have begun comparisons between 
different opacity calculations. 

Two other processes are considered in the stellar 
equations or in stellar pulsations. Both require a de- 
tailed description of the sub-surface layers: the radia- 
tive acceleration and the k mechanism. 



Fig. 2 Influence of the different elements of the iron group 
to the total opacity and comparison between OP and OPAL 
(same representation than Figure 1). 



Radiative accelerations are crucial quantities to 
study the diffusion of elements in stars and to un- 
derstand the stellar photospheric abundances. They 
strongly depend on the atomic properties of the dif- 
ferent ions, so they largely depend on the quality of 
the atomic data and on the knowledge of the ionization 
stage of the different elements. For a detailed descrip- 
tion of this process, we refer to the book dedicated to 



G. Michaud, and in particular to the paper of Alecian 
( 2005 1 . The expression for the radiative acceleration 
is also given in Turck-Chieze et al. (2009). This pro- 



cess must be taken into account in the description of 
stellar envelopes, in addition to the gravitational set- 
tling, turbulence and sometimes magnetic activity. For 
radiative accelerations, one needs to know the detailed 
monochromatic photon flux for all the processes (dif- 
fusion, bound-bound, bound-free...) and for each ele- 
ment. This information is not delivered by the OPAL 



consortium (jiglesias fc Rogers 1995 ) but is available for 



the OP calculations (Seaton & Badnell 2004 1. Some 



specific comparisons lead to 50% differences in the ra- 
diative acceleration for the iron case in the opacity peak 
region mentioned in Figs. 1 and 2 (Delahaye & Pin- 
soneault |[2006| 



It is important to resolve this discrepancy to under- 
stand the chemically peculiar stars and the thousands 
of pulsators observed with COROT and Kepler. The 
rapid variability of the above-mentioned opacities in the 
stellar envelope generates the so-called k mechanism 
and maintains the pulsation in classical Cepheids (due 
to partial ionization of helium) and for /3 Cephei stars 



(due to the iron-group element ionization) (Dupret 



20031 



Calculations for asteroseismology must predict the 
theoretical acoustic modes in order to identify with- 
out ambiguity the observed modes and to understand 
the excitation of the modes. Several papers show the 
diffculties for f3 Cephei stars ( Pamyantnykh| 1999 Bri 



quet et al. 2007 Degroote et al.|[2010[ [DaszynsS^ 



Daszkiewicz & Walczak 2010) and mention that an 



increase of opacity will improve the interpretation of 
the asteroseismic data. 
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Fig. 3 Iron transmission spectra corresponding to temper- 
atures of 15.3, 27.3 and 38.3 eV and densities of 5.5 10"'^ 
g/cm^, 3.4 10~^ g/cm'' and 2.6 10""^ g/cm^ respectively, es- 
timated using the OP calculations. The whole wavelength 
range is useful for the Rosseland mean opacity determi- 
nation (AE — he =12394 A.keV) and corresponds to the 
XUV spectrometer used for the experiment. The ionization 
charge distribution for each condition is superimposed. 

4 Theoretical opacity calculations 



Opacity (cm^/g) 




hnu (eV) 

Fig. 4 Comparison of eight theoretical monochromatic 
spectra for the mean conditions T= 27.3 eV and p = 3.4 
10~'^g/cm^ shown in Figure 3. The hypotheses of these 
computations are explained briefly in the text. 



We investigate the plasma conditions found in the en- 
velopes of /3 Cephei stars for improving the interpreta- 
tion of the asteroseismic observations, without neglect- 
ing previous measurements of iron (and nickel) opacity 
spectra which have successfully led to a revision of opac- 
ities in the nineties, but not always in the range useful 



for astrophysics (Springer et al. 1992; Winhart et al. 
][T995l p^CTTais-Popovics et al. ||2000, ,200l| ). With this 



objective, we perform a comparison between different 
calculations of opacity spectra and an experiment to 
validate these calculations. 

The first step consists in defining a condition com- 
patible with a laboratory experiment. The iron-group 
peak in the envelope of /3 Cephei stars corresponds to 
temperatures of 200,000 K (17.2 eV, 1 eV= 11,600K) 
and densities around 4 10^^ g/cm"^. This density is too 
low to be investigated in a laboratory experiment. So 
we have first looked for conditions which maintain the 
same distribution of ionization stage for iron. 

Figure 3 shows the spectral transmission of photons 
through a heated iron foil of thickness e, T^, = 
exp"'*"'"^, where n^, is the opacity per mass unit in 
cm^/g. Transmission is the quantity measured in the 
experiment. The mean spectrum is obtained using the 
OP tables, for the same distribution of the iron ioniza- 
tion charge (about 9 for T= 320,000 K or 27.3 eV, p 
= 3.4 10~^ g/cm^). Two other spectra also are shown, 
for a variation of temperature of about 40%. In these 
cases, the distribution of charge is displaced by two or 
three units. The figure shows that the spectra vary 
strongly with temperature. In the region where the k 



mechanism is important, the Rosseland mean opacity 
varies rapidly for this reason, so measurement at differ- 
ent temperatures will help to check the whole calcula- 
tion of the pulsation mode excitation. 

Figure 4 shows the opacity spectra versus the photon 
energy obtained by the different calculations generally 
used to describe these plasmas. Clear differences can 
be seen in different parts of the spectral range. The 
corresponding Rosseland mean values show differences 
of up to 40%, in particular for the OP value, but most 
of them agree within 15% despite the very different cal- 
culational approaches. We notice also some variation of 
the mean ionization from one calculation to another of 
about 0.5 units, which is not so small if one considers 
the rapid variation of the opacity with the degree of ion- 
ization, but which cannot explain all of the differences. 
A detailed analysis of these calculations will appear in 
a more detailed paper. We give here a short description 
of the assumptions. 



Based on the average atom model SCAALP (Blan 



card & Faussurier 2004 1, the OPAS code combines 



detailed configuration and level accounting treatments 
to calculate radiative opacity of plasmas in LTE. The 
bound-bound opacity can be calculated using detailed 
line accounting and/or statistical methods. The bound- 
free opacity is evaluated neglecting the configuration 
level splitting. To improve the accuracy of opacities 
into the complex regime of warm dense matter, where 
plasma and many-body effects can be important, the 
free-free opacity is obtained by interpolating between 
the Drude-like opacity and the opacity derived from 
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the Kramers formula including a Gaunt factor and an 
electron degeneracy effect correction. Photon scatter- 
ing by free electrons includes collective effects as well 
as relativistic corrections. 

OP (Opacity Project) is an on-line atomic database 



used in astrophysics (Seaton & Badnell 2004 Seaton 



2007). TOPbase contains the most complete dataset 
of LS-coupling term energies, f- values and photoioniza- 
tion cross sections for the most abundant ions (Z=l- 
26) in the universe. They have been computed in the 
close-coupling approximation by means of the R-matrix 
method with innovative asymptotic techniques. TOP- 
base also contains large datasets of f-values for iron 
ions with configurations 3s'^3p'^3d'', referred to as the 
PLUS-data, computed with the atomic structure code 
SUPERSTRUCTURE. 



The atomic physics SCO code (Blenski,Grimaldi & 



Perrot 1997 Blenski et al. 2000 1 is based on the super 



configuration approach proposed by Bar-Shalom et al. 
[ ( 1989 ) to calculate high-Z element photoabsorption in 
LTE conditions. Here calculations have been performed 
including non-relativistic and relativistic contributions. 

SCO-RCG dPorcherot et al. |2010D is a hybrid opac- 
ity code which combines the statistical super-transition- 
array approach and fine-structure calculations. Criteria 
are used to select transition arrays which are removed 
from the super-configuration statistics, and replaced by 
a detailed line-by-line treatment. The data required for 
the calculation of the detailed transition arrays (Slater, 
spin-orbit and dipolar integrals) are obtained from the 
super-configuration code SCO (see above) , providing in 
this way a consistent description of the plasma screen- 
ing effects on the wave functions. Then, the level ener- 
gies and the lines (position and strength) are calculated 
by the routine RCG of R.D. Cowan's atomic structure 
code ( [Cowan |p8l |. 

The Los Alamos Light Element Detailed Configu- 



ration opacity code, LEDCOP ( |Magee et al.||1995[ ), 
uses a basis set of LS terms based on semi-relativistic 
Hartree-Fock atomic structure calculations (plus a mix- 
ture of LS terms and unresolved transition arrays, some 
of which are treated by random lines for the most com- 
plex ion stages) to calculate opacities for elements with 
Z < 31. Each ion stage is considered in detail and inter- 
actions with the plasma are treated as perturbations. 
Calculations are done in local thermodynamic equilib- 
rium (LTE) and only radiative processes are included. 
In the STA opacity code ( |Bar-Shalom et al. J[T989 ) 



the total transition array of a specific single-electron 
transition, including all possible contributing configu- 
rations, is described by only a small number of super- 
transition-arrays (STAs). Exact analytic expressions 
are given for the first few moments of an STA. The 
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Fig. 5 Transmission spectra of the calculations for condi- 
tions mentioned in Figure 3 and a XUV spectrometer reso- 
lution of 1 eV. Superimposed is the weighted shape of the 
Rosseland mean opacity (dot-dash line). 



method is shown to interpolate smoothly between the 
average-atom (AA) results and the detailed configura- 
tion accounting that underlies the unresolved transition 
array (UTA) method. Each STA is calculated in its own 
optimized potential, and the model achieves rapid con- 
vergence in the number of STAs included. 

CASSANDRA is an AWE (Atomic Weapons Estab- 



lishment, Aldermaston, UK) opacity code (Crowley et 



al. 2001 1 used to model plasma in local thermal equilib- 
rium. CASSANDRA'S self-consistent field calculation 
uses the local density approximation for bound states 
and has a free electron contribution based upon the 
Thomas-Fermi model. The Cassandra data are British 
Crown Owned Copyright 2010/MOD. 

FLYCHK is a free access atomic spectroscopy code 
( Chung et al.||2005 ) used to generate atomic level pop- 
ulations and charge state distributions and synthetic 
spectra for low-Z to mid-Z elements under LTE or 
NLTE conditions 'http://nlte.nist.gov/FLY/, 



5 Experimental challenges 

The new spectral opacity measurements have been per- 
formed on LULI 2000 www.luli.polytechnique.fr/ with 
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two lasers in 2010. A nanosecond laser delivers an en- 
ergy between 30-500 J in a 500 ps duration pulse. This 
laser is used to irradiate a gold cavity (hohlraum) on 
which an absorbing foil of the considered element is 
deposited and heated. After a required delay to get 
the proper density and temperature, a picosecond laser 
interacts with a backlighter to produce x rays which in- 
teract with the formed plasma: for a description of the 



experimental setup, see Loisel et al. (20091. We detect 



the transmitted opacity spectrum on a streak camera 
placed behind a newly designed XUV-ray spectrometer. 
Obtaining high-quality measurements required address- 
ing many challenges based on experience with previ- 
ous measurements. A detailed analysis of the previous 
experiments and their difficulties is found in [Chenais^ 
Popovics ( 2002 1 ; Bailey et al. ( 2009 ) . We summarize 



some of these relevant for astrophysical applications: 

- One needs to form a plasma in LTE. This require- 
ment supposes a good simulation of the whole exper- 
iment to probe the foil in the best conditions. The 
rapid expansion of this foil during heating is limited by 
placing the foil between two thin samples of a low-Z 
material (typically carbon). 

- We measure the foil temperature with a /^Dmix 
and limit the temperature gradient by placing the foil 
between two cavities. The temperatures depend on the 
geometry of the experiment, the position of the foil, and 
the correctness of the simulation. 

- The resolution of the XUV spectrometer is impor- 
tant to discriminate between calculations; it varies with 
the wavelength and one must avoid saturation. 

Figure 5 shows the transmission obtained for the dif- 
ferent calculations including a spectrometer resolution 
of 1 eV for the the conditions mentioned in Figure 3; the 
differences must be distinguishable by the experiment. 
Notice the distinct behavior of the OP calculation at 
low energy and some saturation of the spectra which 
illustrates the difficulties of exploring the low-energy 
regime. 

To reduce the difficulties noted above, we measure 
the spectra of 3 or 4 neighboring elements (chromium, 
iron, nickel, copper) using two thicknesses for each of 
them. Moreover we use two cavities to heat the foil on 
the two faces to limit the temperature gradient. 

In summary, this study shows that it would be of 
interest to develop complete opacity tables including 
photon energy spectra with the best physics available 
for asteroseismology. It is important to improve the 
coverage of some crucial regions for interpreting the 
space observations. Some special validations by lab- 
oratory experiments are required to validate "the best 
physics" . 
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